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Abstract 

A generalized model of Heisenberg quantum antiferromagnet on an 
arbitrary graph is constructed so that the VBS is the unique ground 
state. The norm of the base state and equal time multi point corre- 
lation functions are computed in terms of generalized hyper geomet- 
ric functions. For the one-dimensional periodic Heisenberg model we 
present a method of computing multi point correlation functions based 
on the study of a commuting family of transfer matrices. The connec- 
tion of multi point correlators with Young tableaux and Gegenbauer 
polynomials is found. 

Introduction. 

Theory of antiferromagnetism is very important. We consider the models 
with Valence Bond Solid ground state [2]. The study of generalized Heisen- 
berg antiferromagnets (see [3]-[9]) is of great interest . A distinguishing 
feature of these models is that their Hamiltonians can be represented in the 
form of a linear combination of projections, which makes it possible to ex- 
plicitly construct the base states for the models in question. The models we 
study are a generalization and modification of those considered in [2], [5], [6], 
[7], [8]. They have a valence bound state (abbreviated VBS) as the ground 
state (which is distinct from the Neel ground state). In the first part of the 



work - §§1-7 - we devote our attention mainly to questions of uniqueness of 
the VBS. We prove that the VBS is the unique ground state of the model 
in any dimension for a periodic lattice with any coordination number. Thus, 
the Neel state is never a ground state for the model in question. We have 
succeeded in reformulating the model on an unclosed chain (by introducing 
special boundary conditions) in such a way that the VBS is, as before, the 
unique ground state. We can construct a generalized model of a quantum 
antiferromagnet (and boundary conditions) on an arbitrary graph, so that 
the VBS is the unique ground state. 

The second part of the paper - §§8-13 -is devoted to the computation 
of the square of the norm of the VBS of the wave function and of multi 
point correlators for the 0(m)-Heisenberg model on an arbitrary graph. The 
computations make essential use of the properties of harmonic polynomials 
A"^"(n) (Theorem 7) introduced in §9, which are a natural generalization 
of Gegenbauer polynomials to the case of several variables. Theorem B in 
Appendix B plays a central role; it makes it possible to derive recurrence 
relations for the generating function of the correlators (Theorem 9). We show 
that the Heisenberg 0(m)-model on one-dimensional or periodic chain are 
connected with a commuting family of transfer matrices (§9), which provides 
an alternative method of computing the correlation functions (§§10, 11). 

An outline of our paper is as follows. In the first two sections we present 
a general construction of models having a unique ground state. The models 
describe the intersection of quantum spins-distinct spins are situated at the 
vertices of an arbitrary graph F. Further, for the models constructed we 
explicitly give the Hamiltonian and present the construction of the ground 
state. Generally speaking, the ground state is not unique. In order to for- 
mulate a condition for the VBS, we introduce the required notation. We 
denote by A^o the number of vertices and by A^i the number of edges of the 
graph F. Let Si be the value of the spin situated at the vertex of F with 
index /. We consider the vector 2S whose components with index / is equal 
to 2S. We denote by (kl) the edge of F joining the vertices k and /. With 
each edge (kl) of F we associate a number Mki which we henceforth call an 
alternating number. We now consider the vector M of dimension equal to 
the number of edges A^i of F, whose components with index (kl) is equal to 
Mj^i. The (edge- vertex) incidence matrix / [10] is an important geometric 
characteristic of F. This is a rectangular matrix of dimension A'o x A^i whose 
element with index (n, {kl)) is equal to 6n,k + Sn,i- In the second section we 
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prove the following theorem. 

Solvability in nonnegative integers M^i (for fixed spins si) of the system 
of linear equations 

2S = 7-M (1) 

is a necessary and sufficient condition for uniqueness of the ground state 
of generalized Heisenberg magnet corresponding to spins si and alternating 

numbers M^. 

In §§3, 5-8 we study conditions for solvability of the system (1) for a 
fixed collection of spins si. This problem naturally decomposes into two 
problems: i) solvability of the system in integers (the Mki must be integers); 
ii) solvability of the system in positive integers M^i. We give a complete 
solution of the first problem. The answer depends in an essential way on 
whether the graph T is bipartite [10]. We recall that a bipartite graph is 
characterized by the fact that the set of its vertices can be decomposed into 
two subsets A and B so that any edge in F joins only vertices of different 
subsets. In §7 we prove the following theorem. 

For a connected, bipartite graph T the system of equations (1) is solvable 
in integers Mki if and only if the following relation is satisfied: 

E = E (2) 

aeA beB 

For a nan bipartite graph T in ^8 we prove that the system (1) is solvable 
in integers M^i if and only if the following condition is satisfied: 

E^^eZ. (3) 

To prove the theorems formulated above it is useful to first consider cases 
where the graph F is a tree (§§3, 5) or a cycle (§6). For such graphs we solve 
the system (1) exphcitly. The result obtained are used in the proof of the 
theorems in the general case. 

As concerns solvability of the system (1) in nonnegative integers, we 
present only necessary conditions (§6). The question of sufficiency of these 
conditions remains open. 

We consider the model of an antiferromagnet in a quasicrystal separately 
(§4). The fact of the matter is that the thermodynamic limit for models in a 
quasicrystal is analogous to the thermodynamic limit in crystal. For a one- 
dimensional quasicrystal we explicitly compute the multi point correlation 
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functions. In the second part of our work we study multi point correlation 
functions for tiie Heisenberg model on graph F. We begin by considering 
the Heisenberg model on one-dimensional lattice. In this simplest example 
all spins at the nodes of the lattice are equal to 1. The Hamiltonian of the 
model is [5] 

^ = ^ E {s(0 • s(/ + 1) + ^(S(0 • s(/ + 1))^ + 1]. (4) 

Here the S(/) are the quantum spins (a representation of the algebra 0(3); / 
is the index of a node of the lattice). The ground state of the model - the 
valence bond state - is described in §1. The equal time two-point correlator 
is computed in [2] (in the thermodynamic limit): 

{S'^-'{h)-S^^{k)) = ls:i{-3)-^\^^-^^\. (5) 

Here the S"' are the components of the vector S (a = 1,2,3). We prove 
that for the model under consideration the equal time multi point correlator 
reduces to the product of the two-point correlators (5). More precisely, we 
have the following: 

1) If the number of spins is odd, then 

= 0, N = 1 (mod 2). (6) 




We suppose here and below that all the coordinates Ij are distinct. 

2) If the number of spins N is even, 1^ > In -2 > ■ ■ ■ > h > h, then we 
prove that 

/ N \ N/2 

{Es^'iij)) = Uis^^'-M-s'^^-^ikr-i)). (7) 

\j=l I r=l 

A generalization of the model (4) to the case of higher spin S is given in [5], 
[8]. In this case the Hamiltonian of the model is given by formulas (19) and 
(21) of the text. The equal time two-point correlator (in the thermodynamic 
limit) for the Heisenberg model of higher spin S is computed in [5]: 

(5-(Z05"^(i2)) = ^-^^Z (^) ■ (8) 
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In the present paper we compute the equal time multi point correlator 
(in the thermodynamic limit) for a model of spin S. The precise formulation 
is presented in §11 and differs basically from (7). As an example, we present 
a formula for the four-point correlator 



Here > L, > I2 > h. 

We note that in the case of a one-dimensional, periodic chain the models 
considered (4) and (19)-(21) have a unique ground state. For a finite open 
chine the ground state is no longer unique. In the first part of the present 
paper we analyze conditions under which the models considered have a unique 
ground state. For example, for the model (19)-(21) on an open chain the two 
boundary spins s must be replaced by spins equal to s/2. These modified 
models have a unique ground state, such modification does not affect the 
thermodynamic limit, but it simplifies the computation of the correlation 
functions for finite chains. 

In the first part of the paper we prove that there exists a natural, inho- 
mogeneous generalization of the models considered above (distinct spins are 
situated at the vertices of an arbitrary graph) with presentation of unique- 
ness of the ground state. One-dimensional quasicrystals are examples of 
such models. We also compute the multi point, equal time correlators (in 
the thermodynamic limit) for an inhomogeneous model on a chain (see §11). 
In [2] it was shown that the quantum models considered above are equivalent 
to a one-dimensional, classical, modified Heisenberg model, which is useful 
in computing the correlators. Wc give a description of the modified model in 
§8. Using the integral operators Km (sec §9), we construct a transfer matrix 
for the classical model and find the spectrum and eigenfunction of the oper- 
ators Km in §9. This enables us to compute the norm of the wave function 
(see §10) and also to find the multi point correlators (§11). In the following 
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§12 we give a generalization to the case of the group 0{m) of the modi- 
fied Heisenberg model considered earlier. In §13 we consider the quantum 
and corresponding classical Heisenberg model for a multidimensional lattice. 
The inhomogeneous Heisenberg model (arbitrary spins are situated at the 
vertices) for a complete graph (all the vertices of the graph are joined by 
edges) is a natural algebraic object. We consider the generating function for 
the correlators and derive for it recurrence relations which make it possible in 
principle to find it for the complete graph (and hence for an arbitrary graph). 
We show that the generating function of the correlators can be expressed in 
terms of generalizes hyper geometric functions. 

We consider in more detail a one-dimensional quasicrystal [12]. There 
exists a quasi periodic covering by a direct infinite sequence of two intervals 
(short and long). One of the examples of such a covering is connected with 
the golden section r = (1 + ^/5)/2. The position of the end of the Ith segment 
of the covering is found by the formula 

(10) 

Here a is real parameter. We denote by pl {Ps) the probability of the occur- 
rence of long (short) interval [12]. We construct a Heisenberg antiferromagnet 
for the quasi periodic lattice (10). To each long segment we assign an alter- 
nating number and to each short segment we assign a number Mg. The 
spins (36) can then assume only the three values 

s{l) = {Ml, Ms, {Ml + Ms)/2} (11) 

depending on the position of the vertex I. The Hamiltonian of the model is 
given by formulas (37), (38), while the ground VBS state is given by formula 
(39). In §10 we explicitly compute the correlation functions - for a finite chain 
and in the thermodynamic limit. We present the answer for the asymptotic 
of the two-point correlator of a quasicrystal: 

{S'^Hh)S''%k)) = ls:i{s{h) + 1)(-1)'-'^ e^p{-m{k - h)} ■ (12) 
Here {I2 — h) 00, 



Xi = l + 



+ a 
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We are deeply grateful to L.D. Faddeev and N.Leskova for his constant 
interest in the work and for many useful remarks. The first version of this 
paper was published as a preprint in 1988 [1]. 

§1. A generalized model of an antiferromagnet. Uniqueness of 
the ground state. 

We consider a periodic lattice in £)-dimensional Euclidean space. The 
interacting quantum spins Sj are situated at the nodes of the lattice (/ is the 
index of a node) . At each node the spins are a representation of the algebra 

su{2y. 

[S^,Sl]=t5ike'"'^S^. (14) 

Here the spin index a assumes the values a = 1,2,3. The lower index I 
numbers the nodes of the lattice. We express the spins in terms of two inde- 
pendent, canonical Bose fields on the lattice - a; and bi. Their commutation 
relations are standard: 

h,a+]=4, [bi,bt] = Si, [ai,bk]=0. (15) 

The components of the spin can be expressed as follows: 

S+ = albi, = blai, 

Sf^^iatai-btbi). (16) 
The value of the spin at the node with index I is an eigenvalue of the operator 

Si^^(atbi + btbi). (17) 

(In a quasicrystal different spins are present at different nodes.) 

Thus, the spin operators (16) act in Fock space whose vectors have the 
form 

P(a+,6+)|0). (18) 

Here P is a polynomial in the variables af' and bf' (/ runs through all nodes 
of the lattice). We seek eigenfunction of the operator S;: 

S;P(a+,6+) I 0) =s;P(a+,6+) | 0). 
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This relation means that the polynomial P is a homogeneous function of 
the variables a'l and b^^ (for given /) of degree 2si, i.e., the polynomial P can 
be represented in the form 

2si 

Here the polynomial does not depend on the variables at the /th node. 

The Hamiltonian of the model describes the interaction only of spins 
situated at the nearest nodes of the lattice (we denote this by {kl) ): 

H^^H(k,l). (19) 

(kl) 

The Hamiltonian density H{k, I) is usually described in terms of powers of 
the scalar product of spins (Sfc ■ S/)". Instead of this we use a special basis 
of polynomials in (S^ ■ S/). These arc the projections Ilj{k,l) onto the state 
with fixed spin J. They can be found from the following system of linear 
equations: 



J=\sk-si\ 



^J{J+'i-)-^-'^i{si + '^)-^Sk{sk + l) 



(20) 

71 = 0, 1, 2, ... , 2Siiiin. 

Here Sk and si are the magnitudes of the spins at the nodes k and I, while 
Smin is the Icast of these two values. The system of equations (20) can easily 
be solved for projections Ilj{k, I); we find the following expression: 

"^*'*= n jj^^^", -ir, - (21) 

jW, \sk-si\<j<sk+s, '^K'^ JKJ^^) 

Here = {Sk + SiY. The projection (21) is a polynomial in {Sk, Si) of degree 
2sniin- In order to produce the Hamiltonian of [1], [3], we suppose that the 
magnitude of spin at each node is the same and equal to s (we give up this 
assumption in next section). The Hamiltonian density is 

H{k,l)= Y AjUj{k,l). (22) 

J=2s+1-M 
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Here M is a positive integer, 1 < M < 2s. (We shall see below that this is 
an important parameter of the theory.) It is important that the following 
relation be satisfied: 

2s = zM. (23) 

Here z is the coordination number of the node (the number of nearest neigh- 
bors). The coefficients Aj are positive real numbers-parameters of the model. 
The model has thus been determined, we call it the AKLT model [2]. We 
present examples. For a one-dimensional, periodic lattice z = 2, and for the 
least value s — 1, M = 1 

H^Aj2{{Si- Skt + 3(S/ • Sfe) + 2} . (24) 

(kl) 

In the two-dimensional case on a hexagonal lattice z — 3, and for the least 
value s = 3/2, M = 1 

= {(S. • S.) + ^(S, • S,)^ + ^(S. . S,)^} . (25) 

(kl) 

We now return to the Hamiltonian (22) and construct the ground state. 

Theorem 1 The Hamiltonian (22) has a unique ground state of the follow- 
ing form: 

\^)^constl[{atbr-arbtr. (26) 

(kl) 

Proof. It is obvious that H > and also that Hj(/c, I) > 0. Thus, if there 
exists a solution IV') of the equation 

H\^) = (27) 

then it is ground state of the Hamiltonian (22). Due to positivity equation 
(27) is equivalent to the set of equations 

Uj{k,im^O (28) 

for any pair of nearest nodes (kl) and for any J in the interval 2s + 1 — M < 
J < 2s. This implies that in adding spins Sk and Si there are no projections 
onto states with complete spin J, where 2s + 1 — M < J < 2s. We now use 
the following theorem, whose proof we present in Appendix A. 
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Theorem 2 (on addition of spins). We suppose that after addition of spins 
Si and Sk there arises a state \tlj) with zero projections onto spins J of the 
interval 

si + Sk + l-M <J <si + Sk. (29) 
Then the polynomial is divisible by {a'^bf' — a'^'b'^)^ . 

We continue the proof of Theorem 1. We recaU that we are seeking 
a ground state \'4>) in the form (18). It foUows from Theorem 2 that the 
polynomial P(a+,6+) is divisible by {a'^b'i — a'lb'l)^ for each edge {kl). 
Hence, 

IV') = P{a\ 6+)|0) = \\{{atbt - afbtr] P{a\ b+)\0). (30) 

[{kl) ) 

Here P{a'^, b'^) is another polynomial. We now compute the magnitude of 
the spin of the state 1-0) at the node with index m. This is an eigenvalue of 
the operator 

-^m = 2 («m«m + b^bm) , (31) 

in other words, — Smli^)- 

/,From this it follows that the polynomial P is a homogeneous function 
of the variables a+ and 6+. Applying the operator (31) and (30), we obtain 

2sm ^Mz + Am. (32) 

Here > is the degree of homogeneity of the new polynomial P in the 
variables and b^- Comparing (32) and (23), we find that = 0, i.e., the 
polynomial P does not depend on the variables and b^- Thus, P =const 
as asserted. This means that equation (27) has a unique solution 

m = m4bi-atbtr. (33) 

(kl) 

The theorem on the existence and uniqueness of a ground state of the Hamil- 
tonian (22) has thus been proved. 

■ 

The wave function (33) realized the valence bound state. In [1], [2], [5], 
[6] it is called a VBS (valence bond state). 
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We shall now try to construct a model for an unclosed chain so that 
the uniqueness theorem is preserved. We thus consider a one-dimensional 
chain of N nodes I — 1,2,. . . ,N. We observe immediately that for the 
interior nodes (/ = 2, . . . , — 1) the coordination number z = 2, while for 
the boundary nodes (/ = 1, N) z = 1. From the relation 2si = ZiM (for 
M = 1) it follows that for interior nodes Si = 1 while for the boundary nodes 
sn = 1/2. 

We take the Hamiltonian in the form 

X;n.,+.,+,(M + i)- (34) 

1=1 

We note that for interior edges 11 is the projection onto spin 2: ITsj+sj^^ = 112, 
while for boundary edges it is the projection onto spin 3/2 (11^;+^,^^ = 113/2 
for I — 1 and / = N — 1). This is the difference between our Hamiltonian 
and that proposed in [2], [5], [6], [7]. In those papers it was assumed that 
Si — Sn — 1, and the ground state is fourfold degenerate. For the hamiltonian 
(34) it is easy to prove that the ground state is unique. In the next section 
we construct a model with a unique ground state on an arbitrary graph. 

§2. The model on an arbitrary graph. 

We shall construct a generalized model of an antiferromagnet on an arbi- 
trary graph r. We do this so that the theorem on existence and uniqueness 
of the wave function of ground state (VBS) is preserved. We consider an ar- 
bitrary graph. It consists of A^'o vertices. Some of the vertices are connected 
by edges. We denote the number of edges by A^i. The edge-vertex incidence 
matrix / plays an important role below. This is matrix with A^'o rows and 
Ni columns. Each of its matrix elements is equal to or 1. Wc define / 
more precisely. Each row of / is connected with specific vertex of the graph 
F; each column is connected with an edge. If a vertex belongs to an edge, 
then the corresponding matrix element is equal to 1; otherwise it is equal 
to 0. We have thus defined the incidence matrix. Wc number all vertices of 
the graph, for example, by the letter /. Wc number edges by pair of letters 
(for example, {kl)) denoting the vertices belonging to the given edge. We 
begin the construction of a generalized model of an antiferromagnet. To each 
vertex I we assign a spin si (the integer 2si > 0), and to each edge (kl) we 
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assign another positive integer Mki > 0. At each node these numbers must 
be connected by relation 

2si = Y.Mki. (35) 

{kl) 

Here the summation goes over all edges abutting the node I. This relation 
can be written with the help of the incidence matrix as 

2S = 7-M. (36) 

Here S is an A'^Q-component vector (its components are equal to S/), while M 
is an iVi-component vector (its components are equal to Mki). We discuss 
the solution of equation (36) below, but now we continue the construction 
of the model. We place a quantum spin si at each vertex of the graph. We 
construct the Hamiltonian describing the interaction of nearest neighbors 

{kiy. 

H = J2H{k,l). (37) 

(kl) 

Here the summation goes over the edges. The Hamiltonian density is 

H{k,l)^ ^ Aj{k,l)Uj{k,l). (38) 

J=Sfc+S(+l-Mfc; 

Here Hj is the projection (21) while the coefficients Aj{k, I) are real and 
positive-parameters of the model (depending on J and on the edge). The 
model is thus determined. In analogy to §1 we prove the existence and 
uniqueness of an eigenfunction describing the ground state. This VBS state 
is such that 

where 

(kl) 

We anticipate that Hamiltonian (38) has a gap in the spectrum; the 
correlators (see §11) decay exponentially for quasi periodic coverings, M^i > 
1. 

Several subsequent sections are devoted to the solution of equations (35) 
and (36). The exphcit form of the wave function (39) shows that nuUification 
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of any of the numbers Mu is equivalent to the absence of an edge. It is 
therefore a question of solving system (35) in positive integers M^i > 1. Of 
course, the integers M^/ can be prescribed arbitrarily, and the values of the 
spins can be computed. It is, however, interesting to solve the inverse problem 
and clarify what restrictions the system (36) imposes on the permitted values 
of the spins. For example, it is clear that the spin can be equal to 1/2 only 
when a vertex has one nearest neighbor. We shall consider only connected 
graphs. For disconnected graphs the problem reduces to several independent 
problems. 

§3. Linear graphs. 

We shall first solve system (36) for hnear graphs. We consider a non 
closed, one-dimensional chain with N nodes. Let / be the index of a node, 
I — 1,2, . . . ,N . We complete the relations (35) at the ends as follows: 

Mqi — and M7v,Ar+i = 0. (40) 

A solution of equation (35) has the form 

M,,,+i = i:(-l)'"'-2^fe- (41) 
fe=i 

For I — N it follows from this that 

N 

E(-l)'=-2.. = 0. (42) 

fe=i 

^From the positivity of M it follows that 

M,,,+i = X:(-l)'-'-2sfc>l. (43) 

k=l 

We now consider a periodic chain with an even number of nodes N. The 
equation 

2s, = Mi^i+i + Mi^i^i (44) 
has a non unique solution. A solution of the homogeneous equation 

= M,,,+i + Mi_,^i (45) 
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has the form 

M;,/+i = {-l)'a. (46) 

The nodes with indices 1 and + 1 are identified. This can be used to break 
the chain, i.e., to set, for example, 

Mi,^ = 0. (47) 

This reduces the problem to an unclosed chain, which has already been 
solved. Thus, a solution exists only in the case 

E(-l)' • = 0. (48) 
1=1 

It has the form 

k=l 

Here Mi^n is an arbitrary positive integer. For odd / the left side of (49) is 
positive, which imposes restrictions on the right side of (49). Thus, recalling 
positivity, it is possible to change {M; ;_|_i} for fixed values of the spin s. 
This means that for the same choice of spins it is possible to produce several 
distinct Hamiltonians (38) whose definition contains the collection {Mi^^i}. 
Each of these Hamiltonians will have a unique eigenfunction of the ground 
state. 

We now consider a cycle with an odd number of nodes. Equation (44) 
has the unique solution 

I N 

fe=l k=l 

N 

^From this it is clear that the sum E(~-'-)*^ ' ~ (— l)^-/^i,iv must be an 

k=l 

integer. Since we have the equality 

N N 

Y^Sk + j:i-^)'-^k^2 ^ Sk, (51) 

k=l k=l k=0 (mod 2) 
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N 

it follows that the sum ^ sj. is an integer. Hence, 

k=l 



N 

2sk = an even number. 

k=l 



(52) 



§4. Quasicrystals. 

There are now a large number of works devoted to quasicrystals (see, for 
example, [11]-[14]. Analysis of dynamical system in quasicrystals is also of 
broad interest. For example, it has been possible to solve the Ising model 
and the eight-vertex Baxter model in two-dimensional quasicrystal [14] . 

In connection with the problem of finding the densest packing specific 
associated with quasicrystals are of major interest [15]. 

Here we shall construct a generalized model of a quantum antiferromag- 
net in a quasicrystal of two or higher number of dimensions. The simplest 
example of a quasi periodic covering is the Penrose tiling of the plane by 
translations of ten rhombuses. Corresponding figures can be found in [13], 
[14]. In the situation of general position a quasicrystal (its vertices and edges) 
is bipartite graph [16]. 

An important characteristic of a quasicrystal is the number of nearest 
neighbors of vertices (the coordination number zi). It can run through only 
a finite number of values. For example, for the Penrose tiling of the plane 
by translations of 10 rhombuses we have [15] 3 < < 7. For filling out 
three-dimensional space by translations of 20 parallelepipeds of special form 
(rhomboids) we have [15] 4: < zi < 12. This filling has the symmetry of a 
right icosahedron and reahzes the crystal structure of the rapidly cooled alloy 
AleMn [17], [18]. A generahzed quantum antiferromagnet can be constructed 
in a quasicrystal in the same way as on an arbitrary graph. The general 
outline was described above. We shall present only the simplest example. 

At the vertices of the quasicrystal we place spins which, generally speak- 
ing, are different in magnitude. The spin situated at a given vertex is equal 
to half the correlation number: 



1 



(53) 
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For Penrose rhombuses 3/2 < si < 7/2, while for the icosahedral filUng of 
space by rhomboids 2 < < 6. The interaction Hamiltonian for these spins 
can be taken, for example, in the form 



The Hamiltonian density is the projection onto the highest possible value 
of the spin arising in adding spins at two neighboring nodes k and /. We 
note that the model constructed here of a quantum antiferromagnet in a 
quasicrystal has a unique ground state (the valence bond state) , which differs 
from the Neel state. 

A quasicrystal is a quasi periodic filling of all space. An important char- 
acteristic of it is the following. Any finite part of quasicrystal has an infinite 
set of copies, and these copies repeat throughout space with a particular 
probability. It is just this that makes it possible to justify the presence of 
a thermodynamic limit in quasicrystal. In analogy to [5] it is possible to 
show that there is a gap in the spectrum of the Hamiltonian (54), and the 
correlators of the spins decay exponentially. In analogy to [4] it can be shown 
that the model of quantum antiferromagnet in quasicrystal is equivalent to 
a model of classical statistical physics in the same quasicrystal with Hamil- 
tonian density — log((l — njfcn/)/2). Here n is a unit vector on the sphere. 

§5. Tree graphs. 

For a disconnected graph the system of equations (35) decomposes into 
several independent systems. Therefore we henceforth consider only con- 
nected graphs. The procedure of "cutting-off branches" is useful in the 
investigation of the system (35) for an arbitrary graph; we proceed to a 
description of it. Here an important role is played by bipartite graphs, by 
definition, the set of vertices {1} of bipartite graph can be broken into two 
non intersecting subsets {/} = {A} U {B} so that an edge of the initial graph 
joins only vertices of different subsets. To each vertex / of the graph F it is 
possible to ascribe a parity 



(54) 




(55) 
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We now consider a graph which decomposes into two disconnected subgraphs 
Fi and r2 when one edge is removed, here Fi is a bipartite graph, we denote 
vertices belonging to the edge 6 by a and c (a e Fi, c e F2). It is easy to 
compute the integer M corresponding to the edge b: 

M,= Y.ei- 2si. (56) 

The contributions of the edges of the bipartite graph cancel on the right side 
of (56). The signs are chosen so that €a — I. We now redefine the magnitude 
of the spin at the nodes a and c: 

2sa^2sa-Mb; 2sc^2sc-Mh. (57) 

Thus, solution of equation (36) for the entire graph has reduced to the solu- 
tion of two independent equations of the type (36) for each subgraphs Fi and 
F2 which are not connected with one another. This procedure of reducing the 
system (35) to two simpler independent systems is naturally called the proce- 
dure of "cutting-off branches". We apply this procedure to the investigation 
of system (35) foe the tree graphs. Thus, we consider the graphs without 
cycles (tree graphs). We denote them by the letter V. It is well known that 
any tree is a bipartite graph [10]. We shall study equation (36) in this case. 
We consider any edge (kl). Removing it leads to a decomposition of the tree 
into two disconnected trees T>i U T>2 where k E T>i, I e T>2. 

We choose an alternating sum of spins along the tree "Di, so that Sk enters 
with sign -|-1. It is clear that 

Mki = E • 2sp- (58) 
peVi 

This is a solution of equation (36). it is also possible to express Mki in terms 
of an analogous sum over the second tree. Comparison of the two expressions 
leads to the relation 

Moreover, the following inequality must be satisfied (the condition of posi- 
tivity) : 

Mki^2j2 ^pSp > 1- (60) 
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The Cayley tree is a particular example of a tree. 

We note that in [5] the case where all spins s = 3/2 is considered for a 
finite Cayley tree (with z — 3). It is easy to see that then the system of 
equations (36) has positive integer solutions. The spin must be equal to 3/2 
only for interior nodes, while on the boundary s — 1/2. This guarantees the 
uniqueness theorem. 

We note that it is just the requirement of positivity which does not permit 
taking all spins equal to one another for a tree of general position. 

§6. Bipartite connected graphs. A criterion of positivity. 

We remark that if there is a cycle of even length in an arbitrary graph, 
then a solution of equation (36) is not unique. Indeed, along a cycle it is 
always possible to add to the quantities M a quantity (— 1)'q; analogous to 
(46). It does not change Si. The quantity a can be chosen so that one of the 
quantities M vanished. This corresponds to removing an edge and breaking 
the cycle. Thus, all even cycles can be broken without changing the spins 
but by changing M. For a bipartite graph any cycle is even [10]; therefore, 
the problem on a bipartite graph reduces to the problem on a tree graph. 
Hence, equation (36) is solvable in integers if and only if 

Ee'«' = 0- (61) 

The positivity condition Mki > 1 imposes more complicated restrictions on 
the spin. An obvious consequence of (35) is the condition 

2S; > Zi. 

Here zi is the coordination number of the Ith node. 

For bipartite graphs we shall derive a more refined necessary condition. 
For this we make several definitions. We denote our graph by F. Suppose 
that by cutting n edges it can be broken into two independent subgraphs 
Fi and F2. Suppose that the following conditions is thereby satisfied. All 
vertices belonging simultaneously to Fi and the cut edges belong only to one 
sub lattice, for example, A. We chose an alternating sum of spins along the 
subgraph Fi: 

T.^i-'^^i-T.H Mik. (62) 
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Here ei = 1 for the sub lattice A and e/ = — 1 for the sub lattice B. The sum 
of Mki along the cut edges stands on the right side of (62). 
^Prom (62) it is evident that 

E ■ > n. (63) 

Here n is the number of cut edges. This condition is necessary. We have 
been unable to prove sufficiency of this condition. 

§7. Non bipartite, connected graphs. 

It was shown in §6 that the presence of an even cycle in a graph leads to 
non uniqueness of the system (36). It turns out that two odd cycles joined 
by a chain lead to analogous degeneracy. 




The scheme shown admissible changes of the number M (corresponding 
to edges) which do not change the spins. (Indeed, the sum of the changes at 
each node is equal to 0.) It suffices that the graph shown in the scheme be 
subgraph of F; this already implies degeneracy of the system (36) on F. It 
is easy to see that the graph in the scheme was obtained from an even cycle 
(of length 10) by gluing together two sides. This is a degenerate even cycle. 
Thus, we consider the system (36) on a non bipartite graph. We begin to 
simplify it. We first break all even cycles. We then remove all degenerate 
even cycles so that the connectivity of the graph is preserved. For this we 
remove only non degenerate edges in the operation of breaking cycles in the 
degenerate case. Under this method of breaking degenerate cycles the con- 
nectivity of the graph is preserved [10]. In summary we break all degenerate 
even cycles. After applying the procedure of "cutting-off branches" (see §6) 
we arrive at a connected graph consisting of a single odd cycle. For one odd 
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cycle the problem has already been solved (see (50). The condition for solv- 
ability in integers is the parity of the sum J2 2sa; (see (52)). The condition of 
positivity of M requires further study. 



§8. The modified classical Heisenberg model. 

In [6] it was shown that the quantum Heisenberg model considered in 
§§1,2 is equivalent to a modified classical Heisenberg model. This model can 
be described as follows. The three-components unit vector n(/) (classical) is 
situated at the Ith node of the lattice. The statistical sum of the classical 
model is 

•=/...n-:n:(i-¥^)'" »« 

and is equal to the square of the norm of the VBS wave function 

m-I[i4bt,i-at,ibtr^'^\0). (65) 

l=~L 

The multi point correlation functions are computed by the formula 

N N r ^ 

{^\ us'^i^j) IV') = n(«a.-)+i) L n (66) 

j=l j=l "'■^ l=-L 

' ^ l-n(0.n(i+l) \""'" 



X 



n nn 

l=-L \ ^ / j=l 



Here rjv > r^-i > • • • > T2 > ri, s{l) in the magnitude of the spin at the 
Zth node, and n°'{l) is the component with index a (a = 1, 2, 3) of the vector 
n(0. 

We introduce the integral operator Km- It acts in the space of functions 
/(n) on the unit sphere S'^ according to the formula 

{KMf){n,) = j^jn {^-^^^y f{n,). (67) 

In the next section we diagonalize the operator Km and show that the op- 
erators Km from a commuting family. An individual operator Km is linear 
combination of M-|- 1 projections. Using the properties of the operators Km-, 
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we compute in explicit form the norm of the function (65) and the correlation 
functions. 

We note that the VBS wave function of the inhomogeneous Heisenberg 
model for an arbitrary graph was found in (39), §2. The statistical sum for 
the modified classical Heisenberg model is equal to the norm of the VBS 
wave function (39) and can be computed by the formula 

"''^ I (ki) ^ ^ / 

Here Mki > is an arbitrary choice of alternating numbers. For the complete 
graph all the vertices are joined by edges. If for some edge (kl) we have 
Mki — 0, then this is equivalent to the absence of the edge (kl) in the graph. 
Hence, we can consider an arbitrary graph as special case of a complete 
graph. 

§9. Commuting transfer matrices. 

We consider the integral operator Km with kernel 

KM(ni,n2) = (^^^^). (69) 

We shall first find its eigenfunction. We denote them by A;^(n) - these are 
symmetric, traceless tensors of rank N. For example, 

(70) 

Before writing out the general formula for A^(n), it is useful to introduce 
some notation. For a collection of indices — {aj} we set 

N 

n^^^Yln"^. (71) 
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For even N, N = 2k, an important role is played below by a partition of the 
set An into pairs = ^aidaba}- We define the delta function of the set 
Ar^ by the formula 

= ^ E W:- (72) 

We note that the number of terms in (72) is equal to 

(2A;- 1)!! = (2A;)!/(2*^- A;!). 

Theorem 3 The eigenfunction A^^ of the integral operator Km are given 
by the formula 

[N/2] 

A^-(n) = E ^k{N) E n^^'^^-^SiA,,). (73) 

k=0 AafeCAjv 

Here (Ajv/^2fe) U A2k — A^. The coefficient jk{^) is 

The tensor A"^^ is a polynomial in n" whose leading component is equal to 

Theorem 4 The polynomial A^''' (n) is an eigenfunction of the operator Km, 

Ij^i {^^—^^^YA^-{n^)^z{N,M)A^-{n2), (75) 
with eigenvalue 

1 ^-^ i-M 1 
z(N, M) = TT — ^(0, M) = — — . (76) 

It is interesting to note that 

z(A^,M)=0, if AT > M + 1. 
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Hence, the integral operator Km is a linear combination of M + 1 pro- 
jections. The details of the computation of z{N, M) are presented in Ap- 
pendix B. We note that the A^^(n) are eigenfunction of the Laplace operator 
on the sphere S'^ (for details see §12). 

The of the operator Km do not depend on M (see (73)). This implies 
that the operators Km commute for different values of M: 

[Kmi, KM2] — 0. 

The norm of the VBS wave function or the statistical sum for the mod- 
ified Heisenberg model (64) and the correlation functions (see §10) can be 
expressed in terms of the transfer matrices 

r(r2,ri)= n'xM(o. (77) 

l=ri 

All factors in (77) commute with one another. The set of eigenf unctions of 
the transfer matrix T is given by (73) and does not depend on the collection 
of numbers {M(/)}. Hence, the transfer matrices (77) form a commuting 
family for different collections of the numbers {m(/)}. 

The number of independent components of the tensor A^^ (which is sym- 
metric and traceless) is equal to 2A'" -|- 1. The tensor A^^ generates an 
irreducible representation of the algebra o(3) of spin N. The product of 
irreducible representations decomposes into a direct sum of irreducible rep- 
resentations. We shall need the explicit form of this decomposition in the 
case where one of the spins is equal to 1: 

1 ^ 

Aj^ -n - A^_^^ 2N +1 ^/^"^ "^+1 ^ ^ 

2 ^ ^ 

^ iai...ai...aj...an+icai 

We recall that — {a^, . . . , ajsf}. We note that in adding spin N and spin 
1 only the spins N + 1 and — 1 occur. We further define the projections 
P onto spin A^ -|- 1 and spin A^ — 1 by the formulas 

P^r+i (A^- • n«-+i ) = Ai-^«-+^ , (79) 
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1 ^ 



_ f . ai...ai...aj...ajv+l roi 

(2iV+l)(2iV-l),<4<^ 

We remark that Pat-i + Pn+i = 1 which follows from (79). The value of the 
projections with other indices on the polynomial A^^n"^+i we define to be 
zero: 

Pj(A^^n"^+i) = 0, if J ^N±l. 

It is natural to consider the number to be the spin of the tensor A"^^ . 
In computing multi point correlators it turns out to be useful to consider the 
following combination of projections: 

PoK-Pfe,_,K--iPfc,_,K-- • ■■Pk,{n'''Pk,{n''''PkX'"'))) ■ ■ •)• (81) 

The expression (81) is well defined and is a linear combination of products 
of the delta functions 5"!^. It is different from only if A;2 = 0, 2 and for 
j = 2, . . . , N — 2 the equalities |%+i — kj\ = 1 are satisfied. We present some 
examples: 

Poin"'"') = ^S2l, P2(n"i«^) = A^i"% Po(n"3Pfc(n"i«^)) = 0, 

Pi(A^i"^n"^) = ^{n'^'S^l + n'^'S^l) - ^W^'S^, (82) 

Po(n-Pi(n«3P2(n--))) = ^{S^S^ + S^S^) - ^S^S^. 

A combination of the projections (81) is a collection of nonnegative integers 
{ko, ki, ki, . . . , /cjv-i, kpf} such that kg — 0, \ki+-i — ki\ — 1, i — 0,1, . . . , N — 1, 
k]^ = 0. It is clear that it is also possible to consider such collections with 
another boundary condition = I- It is obvious that I < N, I = N (mod 2). 
In §11 we show that for fixed kN '■= I such sequences can be parameterized 
by standard Young tableaux of the form {{N + l)/2, {N — l)/2). Hence, the 
number of different combinations of projections of the form (52) is equal to 
for A^ odd and equal to the Catalan number Ck — {2k)\/k\{k + l)\ \i N = 2k. 

We finish this section with the remark that in principle, using the rules of 
passing from functions of spin variables to functions of spherical coordinates 
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on the sphere 5*^ described in [19], it is possible to find an expression for A"^^ 
in terms of the original quantum spins S"" (see §1). We illustrate the nature 
of the answer only with one example. We consider Ag^"^ = ^aia^ _ i^ai 
terms of the quantum spins A2 can be written as follows: 

' (s + l)(2s + 3) 

§10. The norm and the two— point correlator. 

We first consider an open chain and compute the norm of the VBS wave 
function (65). We shall proceed from formula (64). We remark that it is 
possible to rewrite the expression for the square of the norm of the wave 
function in terms of the transfer matrix (77): 



^l-n(0-n(/ + l)\''(') 



l=-L l=-L 



L-1 



(83) 



l=-L 



Using (75) and (76), we find 



L-1 J 



The norm of the wave function for a periodic chain can be computed in a 
similar way. Indeed, 

L-1 

{iP I V'),eg = tr n Km(i) = tr T(L, -L). (85) 

I=-L 

All the operators K commute, and their spectrum is known. The transfer 
matrix T(L, —L) has m + 1 nonzero eigenvalues. Here m — min{M(/)}. 
According to (76), these eigenvalues are equal to 

L-1 

Zjv= n Z{N,M{1)), N^O,l,...,m. (86) 

l=-L 
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Hence, the square of the norm of the wave function for a periodic chain is 

m 

I ^)reg= Y.{2N + 1)-Zn. (87) 
Ar=0 

The factor 2A^ + 1 in (87) describes the multiphcity of the degeneracy of the 
eigenvalue Z^- 

An especiaUy simple formula for ('0 | ■0)reg is obtained for m = 1 = M{L): 

{V^ I ^;reg ^^(^^ ^ ^ 2 (M(/) + 1)(M(/) + 2) • 

We now consider the thermodynamic limit for {ijj | ■0)i.eg as L — > oo. It is 
easy to see that 

m L— 1 -1 

(V' I V^)reg = E (2^ + l)^iv - ^0 = n liTTTY-T = I (88) 
iV=o ;=-L -'^^ V'J "I" 

Hence, the thermodynamic limit does not depend on the boundary condi- 
tions. Computation of the limit in (88) is based on the inequality 

Z{N2, M) < Z{Ni, M) if A^i < N2. (89) 

We now proceed to the computation of the two-point correlator for an 
open chain. We use the representation (66): 

^0-7 n d^i n [l:ii^W-^<l±l) ] ""^'^ n- (rO . n- (r,) . (90) 



l=-L 



Here r2 > ri. Formula (90) for the correlator can be written in terms of the 
transfer matrix (77) in the following manner: 

i/jo' ■ Ao • T(L,r2)n"^(r2)T(r2,ri)n'^i(ri)T(ri, -L)Ao. (91) 

In (91) the transfer matrices are defined in analogy to (77). For example, 

L-l 

r(L,r2)= n^M(0, r2<L. (92) 

l=r2 
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The function Aq is an eigenfunction for tlie transfer matrix T(ri, —L) with 
eigenvalue 

It is easy to see that (93) is contained as a factor in {t/j \ ■0)o and hence 
cancels in (91). The function n"i(ri) = A"^ is an eigenfunction for (92) with 
eigenvalue 

-M{1) 

11 (M(0 + l)(M(/)+2)- ^^^^ 

Further, it is clear that T{L, r2) does not contribute to the correlator. Finally, 
for (90) we find 

{n"(r.)»"(n)> = i^Sn(i|^). PS) 

Formula (95) gives an expression for two-point correlator for finite open 
chain. Using (66), we find the correlator of two spins 

Formula (96) gives an expression for the two-point correlator for the quan- 
tum model (38). In the thermodynamic limit formula (86) is preserved (the 
thermodynamic limit does not depend on the boundary conditions). It also 
gives an answer for the one-dimensional quasicrystal (10) in which M(/) 
assumes only the two values Mg and M; in a quasi periodic manner. The 
asymptotic (12) follow directly from (96). If all the numbers M{1) are equal, 
then M = s and (96) reproduces the result (8). 

We now rewrite formula (97) for the square of the norm of the VBS wave 
function in terms of generalized hyper geometric functions pFq. We recall 
their definitions (see [20]): 

J2 («l)fc ■ ■ ■ (ftp)fca^^ ^g^^ 



cti, . . . , Q; 



p 



X 



The symbol {a)k is defined as follows: 

(a)k^^^%-^ ^a{a + l)---{a + k-l), k>0. (98) 
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In terms of hypergeometric series, formula (87) takes the form 

I ^)reg = 




2L+3-r2L+2 



3/2,l,-M(-L),...,-M(-L) 



l/2,M(-L) + 2,...,M(L) + 2 

(99) 

Using Dougall's formula [20] for the completely balanced series 5-F4 we obtain 
an especially simple expression for the square of the norm of the wave function 
for a periodic chain with three nodes (L = 1): 

(V' I V')reg = 

r(M(-i) + i)r(M(o) + i)r(M(i) + i)r(M(-i) + m(o) + m(i) + 2) 
r(M(-i) + M(o) + 2)r(M(-i) + M(i) + 2)r(M(o) + M(i) + 2) ■ 

(100) 



§11. The multi point correlator. 



In §8 it was shown that the correlator for the quantum Heisenberg model 
and the correlator for the modified classical Heisenberg model are connected 
by the relation 

/ N \ N IN \ 

mS-^{r,)) =E(^(0) + l)(n^"^(^.))- (101) 

Here rjy > Tat-i > ■ ■ ■ > r2 > Ti. It is clear that the correlator is equal to 
for odd A^. The right side of formula (101) can be computed by means of 
commuting family of transfer matrices in analogy to the computation of the 
two-point correlator on the basis of (91). The transfer matrix between two 
nearest nodes n"-^{rj) and n"^+i(rj+i) can be by means of the polynomials 

Af.. ^ (we recall that we call kj the spin of the tensor A ''J). The difference 
of two nearest spins kj is equal to fcj+i — kj = ±1 (see (79)). Thus, the spins 
of the tensors form a sequence kg, ki, . . . , oi nonnegative integers such 
that 

ko = 0, \kj+i - %| = 1, j = 0, . . . , AT - 1, kN = 0. (102) 

We use the sequences (102) to write out a formula for the multipoint cor- 
relator in explicit form. Since the function A^^._^ is multiplied by n"^ (rj) at 



28 



the vertex rj, we must use the projections (79). We denote by X[k,M) the 
normahzed eigenvalue of (76): 

^(^'^)=z(m^=(mT^- 

The symbol {a)k is given by formula (98). 

We first write out the answer for correlation functions in the one dimen- 
sional case (all the M{1) are the same, M{1) = M = s; see (22)). 

Theorem 5 We have the equality 

IN \ N-l 

iXln^'ir,)) = Y.X{{X{k,,M)r^^-^^P,{n'^- ...P,^{n^^^^ (104) 

\i=i / {k} 3=1 

The summation in (104) goes over all possible sequences (102). We remark 
that the last factor on the right side of (104) is a c-number equal to a linear 
combination of products of the delta functions 5ai (and does not depend on 
the component of the vector n and the indices of the nodes Vj). It is easy 
to see that 

X(l,M) = ^, X(2,M)= ^(^-^' 



M + 2' ' ' ' (M + 2)(M + 3) 

Using equalities (82) for the projections, wc obtain formulas (8) and (9) 
[s = M). We note that (104) for even N gives an answer for multipoint 
correlation functions also for an open chain (in the thermodynamic limit 
everything remains unchanged). 

We now present an answer for correlation functions in the inhomogeneous 
case (all the M(/) are distinct; see §2). 

Theorem 6 We have the equality 

IN \ iV-l'-j+i-l 

(nn»^(r,))=En n m-,^(own"----n3(^"^^fc.(^""^))---)- 

\i=l / {A;} 3=1 l=rj 

(105) 

The summation in (105) goes over all possible sequences (102). 

The formula for the correlators (12) for quasicrystals follows from (105). 
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We make some remarks regarding the sequences (102) and the projections 
(81). First of all, it is natural to consider sequences (102) of length + 1 
with the boundary condition — I ior fixed I. Let o/^jv be the number of 
them. It is clear that I = N (mod 2). Using induction to N, it is possible to 
show that 

'"-;'^^)-(<"-r-'). (106) 

The full number of sequences (102) of length A'"+ 1 (/cjv is not fixed) is equal 
to 

= ( m ) ■ ^""^ 

^From (106) it is evident that a/ a? is equal to the number of standard Young 
tableaux of the form {{N+l)/2, {N —l)/2) and is also equal to the multiplicity 
of degeneracy of the irreducible representation of the Lie algebra o(3) of sig- 
nature (A'"-|-Z)/2 in the A'"th tensor power of the fundamental representation. 
In our special case I — 0, and hence 

r if = 1 (mod 2) 

ao,N = m\ .^^^ ^ .^^ (108) 
[ A;!(A; + 1)! 

We shall construct a bijection between the set of sequences {k^, . . . ,kN) such 
that /co = ^1 — I, kj e Z+, \ki — /cj+ij = 1, i = 0, 1, . . . , iV — 1, and the 
set of standard Young tableaux T of the form A = {{N + l)/2, {N — l)/2). 
To this end we consider the sequence ko, ki, . . . , k^- We write the number 
i G [1, N] in the first row of the diagram A if — ki = —1 and in the second 
row otherwise. The bijection has been constructed. We present an example. 
Suppose {ki} = (0, 1, 2, 3, 2, 1, 0, 1, 0, 1). Then 

1 2 3 7 9 

4 5 6 8 

We have thus obtained a representation for the iV-point correlator (105) 
in the form of a sum over standard Young tableaux of the form (A^/2, N/2). 
It would be desirable to have a more explicit formula for projection (81), pro- 
ceeding directly from the standard table T of the form ((N + l)/2, 
{N — l)/2). This problem is in the stage of solution. We only make sev- 
eral remarks. We denote by Pt the projection (81) corresponding to the 
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table T. It is clear that Pt is a linear combination of the polynomials A^. 
For example, if 

1 3 ... N-+1 N-l + 1 ... N 

2 4 ... N-l 



then 



if 



then 



. 1 ^ {N-l)/2 ^ ^ (iV-l)/2 

(2AT2) Af-'*'--'" n «r'. 



3 4' 



= (2A + 2)(2A + 4) {'X + KX - \, (110) 



if 



then 



T- 1 2 4 
3 

For the Lie algebra o(3) in the preceding formulas it is necessary to set 
A = 1/2. 

§12. The o(m)— generalization of the modified classical Heisen- 
berg model. 

In the preceding sections we considered the o(3)~modified classical Heisen- 
berg model. There is a natural generahzation to the case of the Lie algebra 
o(m) and irreducible representations corresponding to one-sided Young dia- 
grams. The statistical sum for one-dimensional, classical Heisenberg model 
has the form 

•=t-.n«n(i-^)-". ..... 
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In (111) the vector n G S*™"^, and the invariant measure dQ is presented 
in Appendix B. We emphasize that now n is an m-component unit vector. 
The integral operator Km is defined in analogy to (69): 



M 



{KMf){n2) = dn, ( /(m). (112) 

The eigenfunctions of the operator Km are also given by formula (73) in 
which the coefficient jk{N) are 

The tensor A"^^ (see (73), (113)) is characterized uniquely by the conditions 
of symmetry and tracelessness, and its principal part n"^^ . We note that the 
dimension of the tensor A"^^ (for fixed N) is 



N + m-l\ f N + m-3\ 2N + m-2fN + m-3^ 

(114) 



m — 1 J I m — 1 J m — 3 \ m — 3 

fN+m-2\ fN+m-3\ 



\ m-2 J \ m-2 )' 

and coincides with the dimension of the irreducible representation of the Lie 
algebra o{m) with highest weight (A^, 0, . . . , 0). We shall formulate the basic 
properties of the polynomials A"^^ as a theorem (in which \ — m/2 — 1). 
We consider the polynomials 

21 

m 

Here n e R™, = ^(n")l Hence, A^^ = k^^\s^-i. 

a=l 

Theorem 7 1. The tensor A^'^ is traceless; in other words, if = 
{oi, 02, . . .}; then 



m m 



tr {a„a.}A^^ := E E A^"(n)<^a1 = 0. (116) 

ai=l a2=l 
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2. Let X be a symmetric polynomial in the variables n^,...,n"^, and 
suppose Il{ai,a2}^ = 0- suppose that X'"^^"" = ^ct^^ ■ n"^^. Then 

3. AA"^^ = 0; where A is the Laplace operator for R"*. Hence, A"^^ 
is a harmonic polynomial of degree N , and the polynomials A"^^, Aj^ C 
[1, . . . , m]^ with the condition 4^{ai G : = 1} < 1, form a basis in the 
space of harmonic polynomials of degree N in R™. 

4- Let As be the Laplace operator on the sphere S'^~^ . Then Ash.'^'^ — 
-N{N + m-2). 

5. We consider the operators 

The generators Ejk generate the Lie algebra SO{m): 

[Ejk, Epq] = 5kpEjq, a j <k, p< q. (118) 

The action of the generators Ejk on the polynomials A^'^ (for fixed N) is 
given by the formula 

EjkA^^ = ^ {A(^^^W)\W5a _ A(^^U{i})\{a}^a|_ ^^^9^ 

6. For each polynomial /(n) of degree N, n & R"*, we define its harmonic 
projection (see [19]; here A = m/2 — 1): 

[N/2] j.^i(/^if)(n) 

Then the polynomial A^^ is the harmonic projection of the monomial n^^ : 

A^^(n) = 
Moreover, it is easy to see that 

iSkn^N^2^-k\ (121) 
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7. Let Ai^ = {a, . . . , a}. Then 



N 



Here the are the Gegenbauer polynomials whose definition is presented in 
Appendix B. 

It is interesting to note that formula (73) for A^'^ can be inverted, and 
the monomials n^^ (see (71) can be expressed in terms of the harmonic 
polynomials A"^^: 



[iV/2] 

A 

n 

1=0 A2iCAn 



an^y: r^'m) E ^^""^^'^KMi). (123) 



Here (3i{N) = 7i(M - Z + 1), where 

= 2.(-A-M+lV (^2^' 

The identity (123) plays an important role in further considerations, since 
it makes it possible to solve the problem of decompositing the tensor product 
of two irreducible (single-sided) representations of the lie algebra o{m) into 
irreducible components. In particular, 

1 ^ 

-n - 2(A + N) ^ -^"^ "^+1 ^ ' 

We thus consider projections analogous to (79) and (80): 

Pn-1 ( A^^ • n'^^+i ) = A^^ • - A^^^""^' , (127) 
Pjfc(A^^ • n"^+i) = 0, ifA;^A^±l. 
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We note that the operators Km-, as in the case m = 3, form a commuting 
family: 

[^Mi,^M.] = 0. (128) 

It follows from Theorem 111 that the polynomials A'^'^ are eigenfunctions of 
the operators Km'- 

Ku^ir = ^m(A^, M)Kjf. (129) 
The eigenvalue Z^{N, M) can be computed from the formula 

Z iNM)- r(2A + l) r(M + A + i)(-M)^ 
Z^[N,M)- r(M + 2A + l)(M + 2A + %- 

We have used here the notation (98). The proof of formulas (129) and (130) 
is presented in Appendix B. Another important characteristic is 

which can be computed by means of (130) and is equal to 

Xr^{k, M) = ■ (132) 

^ ' ^ (M + 2A + l)fc ^ ' 

We can now conclude the computation of the multipoint correlation func- 
tions for the one-dimensional, modified, classical Heisenberg (9(m)-model. 
The answer is given by formulas (104) and 105) in which the X{kj,M) are 
now replaced by (132), the projections (79) and (80) are replaced by the pro- 
jections (126) and (127), and the summation both in (104) and (105) goes 
over sequences (/cq, ki, . . . , kj^) such that ko = kn = 0, ki G Z_|_, \ki — ki^i \ = 1, 
i = 0,l, ...,Ar-l. 

In conclusion we present formulas for square of the norm of the VBS wave 
function of the classical Heisenberg 0(m)-model for an open and periodic 
one-dimensional chain (we recall that A = m/2 — 1): 

nin f r(2A + l) ]^'' 'f,' T(M(l) + \ + i) 



We note that if all M{1) e Z+, then 

.-1 

rr 

(2A + 1)m(o 



ii, I v-)„ = n ■ (134) 
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Formulas (133), (134) for A = 1/2 go over into (84): 

(^U) - n r(2A + i)r(M(0 + A + |) 

^^'"^^^^^ " J_i^r(A + i)r(M(0 + 2A + i) ^^^^^ 

/2A,A + 1,-M(-L),...,-M(L) \ 
'■^+^ A,M(-L)+2A + 1,...,M(L) + 2A + 1 i' 



§13. An algebraic approach to the computation of correlators 
for the multidimensional Heisenberg model. 

We consider the statistical sum for the classical Heisenberg 0(m)-model 
on a multidimensional lattice 

4w=L.n<'nn(i^)"- (136) 

sites (ij) 

It is clear that $(M) is equal to the square of the norm of the VBS wave func- 
tion (39) for quantum Heisenberg model describing the interaction of spins 
s — Mz/2 (here z is the coordination number of the lattice). Computation 
of correlation functions of the quantities (1 — n^. • n;)/2, where k, I are nearest 
neighbors, is of major interest. First of all the mean value of (1 — Ui • 10.2) /2 
is the Neel order parameter 

1 - mn^ X _ IUjdnUi^J){{l - • n,.)/2)^((l - n^ • n2)/2) 

Equality to one of the mean of (1 — ni ■ n2)/2 is equivalent to the existence 
of a Neel order in the system. For one-dimensional lattice of 0(m)-spins s 
the Neel order parameter is computed in (134): 



s + 2A + 1 

Here A = m/2 — 1. Equality (138) shows the absence of Neel order for 
the model considered. It is interesting to note that for one-dimensional 
Heisenberg chain there is no correlation between Neel order parameters: 

((1 - n(l) ■ n(2))(l - n(j) ■ n(j + 1))) = (1 - n(l) ■ n(2))(l - n(j) ■ n(j + 1)). 
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To compute the correlation functions of several quantities {l — bfni^-ni)/2 
we consider the generating function for the correlators 

M{M.,}) = / n n (^^^r^ ■ (139) 

sites (ij) 

In the general case all the alternating numbers Mij are assumed to be distinct. 
It is most natural to consider the Heisenberg model for the complete graph 
(all the vertices are joined by edges) with p + 1 vertices 



f /I - n- • n \^'^' 

$.+i({M,})=/ n i-^^) ■ (140) 

j=l i<i<j<p+i V ^ / 



We emphasize again that we consider the alternating numbers Alij as in- 
dependent parameters. In this section we shall prove that $r is a rational 
function of (integer) parameters Mij. We note that an arbitrary graph can 
be obtained from complete graph by annihilating some of the parameters 
Mij. Indeed, the condition Mij = is equivalent to the absence of the edge 
(ij) in the graph. 

We proceed to the investigation of the statistical sum (140) for the Heisen- 
berg 0(m)-model on a complete graph with p + 1 vertices. We wish to find 
recurrence relations between the functions $. For this we first compute the 
integral with respect to a fixed invariant measure dQj of the integral (140) in 
terms of generalized hypergeometric functions. We must thus compute the 
integral 

__,s,,n(i^)-«. (141) 

Before presenting an explicit expression for the integral (141), we give some 
definitions and examples. We consider collections of complex numbers Zij 
and nonnegative integers kij {1 < i < j < p) . We denote these collections by 
z^^ and k^^\ For the collection kS^^ we define 

l<i<j<p l<i<i<P 

(142) 

l<l<i i<l<p 
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We now define the generalized hypergeometric series 



rp{p) I Cili ■ ■ ■ 1 



depending on p(p— l)/2 complex variables z^^ and on p + 1 real parameters 
CKp, /3. The function F^^ is given by means of \p{p — l)-fold sum over 
the collections A;^^: 



[ /3 



Here we use the notation (98): 



l<i<j;<p 



T{a + k) 

r(«) ■ 



(143) 



We note that if the aj are nonnegative integers, then the series (143) 
breaks off and the function 



Trip) I Oil) • • • ) Q^p 



M 



can be represented as a finite sum. We present examples of the functions 
pip) for p = 2,3 (F(i) = 1): 



(2) / ai,tt2 



V- (Qii)fc(Qi2)fc ^ 



(144) 



ir(3) 



«!, q;2, q;3 
/3 



^1, ^2, ^3 



E 



(«i) fci+fc3('^3)A;i+fe2 



A;i+fe2+fc3 



X 



7^2 
^2 



^3 



A;i! k2\ k^V 

We now formulate the result of evaluating the integral (141). 



(145) 
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Theorem 8 Let no, rii, . . . , G S™' ^. Then 

P /l-no-nj\^o^ _ r(2A + l) r(Moi + ■ ■ ■ + Mop + A + i) 



^^(p) / -Moi,-Mo2,...,-Mop 

i -M01-M02 Mo„-A + i 



• (146) 

i<i<j<pj 



^ / r(A+i) r(Moi + --- + Mop + 2A+l) 

zoi — -(W02 — • • • — jwop — A -r 2 I 2 

We recall that A = m/2 — 1. We again emphasize that if all the Mqj e Z+, 
then the series on the right side of (146) contains only a finite number of 
terms. 

^From Theorem 8 we obtain a recurrence relation for the statistical sum 
(140) of the classical Heisenberg o(m)-model on a complete graph. We denote 
the vertices of the complete graph with p + 1 vertices hy I — 0,1, . . . ,p. We 
decompose the set of alternating numbers {Mim} = {Moj} U {Mij} where 
^ ^ h j ^ P- We consider the complete graph with p vertices as a subgraph 
of the complete graph with p + 1 vertices for which the vertex / = and the 
edges (0, 1), ... , (0,p) have been removed. 

Theorem 9 The statistical sums $p and are connected by the relation 



%+i{{Moj}U{Mij}) 



r(2A+ 1) r(Moi + ■ ■ ■ + Mop + A + I) 
r(A + I) ' r(Moi + • • • + Mop + 2A + 1) 



X 



n(-Mo,),(rt 

E —P %{{M,, + k,,}). (147) 

^^^^ (fc(^))!(-EA^o,-A+-)|,(,)| 

The proof follows immediately from Theorem 8 and formula (143). As 
an example we present the formulas for $p for p < 3: 

' ^ ^ r(A + |) r(M + 2A + i)' ^ ' 



X 



n 



[r(A + i)P 
r(M, + \ + \] 



\ r(Mi + M2 + Ms - M, + 2A + 1) ■ 
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It should be noted that Theorem 8 can be used not only for complete graphs. 
We consider, for example, a graph F containing a vertex / of multiplicity 2. 
It may be assumed that I — 1 and the vertex 1 is joined only with vertices 2 
and 3. We denote by F the graph obtained from F by removing vertex 1 and 
the edges (12) and (13). The statistical sums corresponding to the graphs F 
and F are connected by relation 



$r(Mi2,Mi3,M23,...) 



F(2A+1) F(Mi2 + Mi3 + A + i; 



F(A+|) F(Mi2 + Mi3 + 2A + 1) 



(149) 



X 



(-Mi2)fc(-Mi3)fc 

-Mi2 - Mi3 - A + |)fc 



$p(M23 + A;,...). 



The statistical sum $r was computed in §§10, 12 for a periodic, one-dimen- 
sional chain (a polygon). If (135) is substituted into (149) we obtain a new 
identity for generalized hypergeometric series: 



2A, A + 1, —mi, . . . , —nip 

A, mi + 2A + 1, . . . , mp + 2A + 1 



(150) 



X 



X 



F(A + I) F(mp-i + 2A + l)F(mp + 2A + l)F(mp-i + m^ + A + |) 
F(2A + 1) ' F(mp_i + A + |)F(mp + A + |)F(mp_i + m^ + 2A + 1) 

(-mp_i)fc(-mp)fe(A + |)fc 



k>Q 



m. 



.p-\+\\{2\ + l)k 



p+iFp 



2A, A + 1, —mi, .... — mp_2! —k 

A, mi + 2A + 1, . . . , mp_2 + 2A + 1, A; + 2A + 1 



In particular, for p = 3 we arrive at Doudall's formula [20] for the completely 
balanced series 5F4. For p = 4 we arrive at the formula for the transformation 
of the completely balanced series 6-P5 into the Saalschtz series 4F3: 
2A, A + 1, —mi, — m2, — m3, —1714 
A, mi + 2A + 1, . . . , m4 + 2A + 1 

F(A + |)F(mi + 2A + l)F(m2 + 2A + l)F(mi + m2 + A + |) 
F(2A + l)F(mi + A + |)F(m2 + A + i)F(mi + "^2 + A + 1) 



6^5 



X 



—mi, 1712, ^3 + ^4 + 2A 

-mi - m2 - A + |, m3 ^ 



hl,A + i 

2A+l,m4 + 2A + l 



(151) 
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We note that the recurrence relation (149) makes it possible to compute the 
statistical sum $r for the following graph F: 



mi 



r = < 



m4 



m2 



m3 



^ _ [r(2A + 1)]^ r(mi + ma + A + |)r(m3 + A + |)r(m4 + A + |) 
[r(A+i)]4 r(mi + m2 + 2A + l)r(m3 + m4 + 2A + 1) 

r(m5 + A + |)r(m3 + m4 + m5 + 2A + 1) 



X 



r(m3 + m5 + 2A + l)r(m4 + mg + 2A + 1) 



X F ( "^2' "^3 + "T.4 + + 2A + 1, m5 + A + I 

^ ^ y -mi - m2 - A + |, m3 + ms + 2A + 1, m4 + ms + 2A + 1 

The proof of Theorem 8 makes use of a large number of different identities 
for sums of products of binomial coefficients; they are rather cumbersome and 
will be presented in a separate paper. Here we note only that Theorem 8 
follows from the more general Theorem B, whose formulation is presented in 
Appendix B. 

Appendix A. 

We represent the operator a'^ as the operator of multiplication by Xi and 
the operator b^^ as the operator of multiplication by yi. The spin operators 
at the Ith node can then be represented in the form 

oxi dyi 

We consider the space of the irreducible representation of spin s. A basis 
can be taken in the form 

Vu. = xf'^^yl''^; II = -si, si + l,...,si-l, si. {A.2) 
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The spin operators act on the basis in the following manner: 

S^Vi_, = (si - ii)v^+i; 

Sfv,^ = {si + Ai)v^-i; (^.3) 

We now consider the tensor product of two irreducible representations Vg^ <S> 
Vsi^. The spin operators add: 

S^^Sf + Sl S^^S^ + Sl (A4) 

The tensor product decomposes into a direct sum of irreducible representa- 
tions: ^ 

J=\si+Sk\ 

We shall construct the leading vector vj of the irreducible representation Vj. 
We seek it in the form 

— 2^ ^iJ-hiJ^k^i Vi Vk ■ K^-'^) 

IJ'j+IJ.k=J 

The leading vector is a polynomial satisfying the two equations 

S'vj^Jvj- S+vj^Q. (A7) 

These equations reduce to recurrence relations for the coefficients am^nf. which 
have the unique solution 

vj = xf^-^xl"{yiXk - y^xiY. (A.8) 

Here /i — si + Sk — J- 

The remaining vectors of the representation Vj can be obtained from the 

leading vector by means of the action of the lowering operator S~ . How- 
ever, it should be noted that the operator ^-commutes with the operator of 
multiplication by bracket in (A. 8): 

[S^,yiXk-ykXi\^0. {A.9) 
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/,From this it is follows that all the vectors of the irreducible representation 
Vj are divisible by {yiXk - UkXiY- 

We have thus proved that if it is known that after addition of two spins 
si and Sjfc there is no projection onto states with full spin 

Si + Sk+l- H<Si + Sk, 

then all the vectors of the representation are divisible by {yiXk — yk^i)^- We 
recall that a'^ — x^, b'l — yi. 

Appendix B. 

To compute the eigenvalues of the integral operator Km for the Heisen- 
berg 0(m)-model on a one-dimensional chain we use properties of the Gegen- 
bauer polynomials Cj^{x). We recall [19] that 

We assume below that X = m/2 — 1. The relation between the Gegenbauer 
polynomials Cj^{x) and harmonic polynomials A^^(n) defined by formulas 
(73) and (113) consists in the following equahty: 

Y^A^-{n,)nt = :p^CUni ■ n,). {B.2) 
Formula (B.2) can be verified directly. We note that 



Hence, 



We now proceed to evaluate the integrals over the (m — l)-dimensional sphere 
of functions containing the Gegenbauer polynomials C^(nj ■ rij). We first re- 
call (see [19]) the definition of spherical coordinates in R'" and the invariant 
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measure on the sphere S*™ ^. Cartesian and spherical coordinates are con- 
nected by the transformation 

xi — r sin 9m-i ■ ■ • sin 02 sin 9i, 

X2 — r sin 9m-i ■ ■ ■ sin 62 cos 9i , 

iB.5) 

Xm-i = r sin 9^-1 cos 9m-2, 

Xm = rcos^^_i. 

Here < r < 00, < ^1 < 27r, < ^fc < tt, A; ^ 1. 

The invariant measure on the sphere 5*"*"^ is defined by the formula 

= ^—1 • • • ^2C^^i • • • d9m-i. {B.6) 

Using the transition to spherical coordinates, it is not hard to show that 



r(2A + i)r(M + A + l)r{N + a + 1: 

[T{X+WT{M + N + 2X + 1) 



2F1 



-N,M + X+l 



{B.7) 

We recall that A = m/2 — 1. Further, using (B.7) and (B.l), it is possible to 
evaluate the integral 



/ 



l-Ui- n. 



M 



r(2A + l)r(M + A + i)(-M)^ 



{B.8) 



r(A + |)r(M + A^ + 2A + l) 

We now proceed to the proof of Theorem 70 and formula (129). To this end 
we consider the integral 



e(nO 



dQ 



1 — n, ■ n~ 



M 



A^-(n,). 



(S.9) 
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/,From the Theorem 7, 2) it follows that ©(rij) = cA'^^(nj). In order to find 
the constant c, we consider the convolution (B.2) of the function 0(nj) with 
nf^. By (B.8) we find 

_ r(2A+l)r(M + A + |)(-M)^ 

r(A+i)r(M + 7V + 2A + l) ' ^ ' 

which coincides with the eigenvalues Zm{N,M) given by formula (130). We 
formulate the result so obtained: 



Sm-1 V 2 



r(2A + l)r(M + A + |)(-M)^ 



~ r(A + |)r(M + 7v + 2A + i) ^ ^^-^^^ 

We shall say a few words regarding the proof of Theorem 8. The details will 
be published in a separate paper. From identity (B.ll) and the inversion 
formula (123) we obtain the equality 



Sm-l V 2 



'j 



r(A + l)r(M + A+l) ^ ^ ^ 



A21CA 

where, by definition, 



^Prom identity (B.12) it is easy to derive the equahty 



r(2A + i)r(M + A + i) ^ X / 

r(A+i)r(M + A+l)\ i: Zi;'WM)(„..„.„,.J*(n.J, (B.14) 
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Here we have used the notation 



(no • ua) = n(^o • Ha), 



aeA 



SM^l E n(na.-n,J (R15) 

(see definition (72)). We note the following consequence of formula (B.14): 

f f ^ oiv/2 r(2A + l)r(M+f + A + i) 

/ asio(no • n^.,) = 2 ' — r- — ^ ^ d(n4„). (B.16) 

Jsm-i °^ ° r(A+ |)r(M + Ar + 2A + 1) ^ ^ ^ 

The next step is proving Theorem 8 consists in considering the integral op- 
erator Km,.m,- C^iS""-^) C^iS"^-^ X 5'"-!): 

(l?„.,«./)(n. . „.) = da (i^ipi)"' {^^^f /(^). 

(B.n) 

It is clear that the case p = 2 of Theorem 8 corresponds to the evaluation 
of the function {Kmi,M2^){'^i " 112). This function can be most simply com- 
puted by using the expansion of the series ((1 — a;)/2)^ in the Gegenbauer 
polynomials C^(x) for A = m/2 — 1: 

fl-x\^ r(2A + l)r(M + A+l) (iV + A)(-M)^ , 

V 2 J ~ r(A + i)r(M + 2A + l)^oA(M + 2A + %'^^^''^ ^''■'''^ 

(see, for example, [19]). As a result, we obtain 



r(2A + i)r(Mi + M2 + A + ^) 
r(A + |)r(Mi + M2 + 2A + 1) 



2F, 



-Mi,-M2 
-Ml - M2 - A 



ni • ns 



(5.19) 

It is clear that equality (B.19) corresponds to the case p = 2 of the Theorem 8. 

The next step consists in computing the function Kmi,M2{''^o-''^An) ^'^^ 

An. This computation is based on the expansion (B.18) and identity (B.14). 
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We shall not present the formula for this function but rather immediately 
formulate a general result having Theorem 8 as a special case. 

Theorem B. Let no, rii, . . . , e 5™"-^. Then 



1 - no • ni 



Ml 



1 - no • n. 



Mr, 



(no-HA^v) 



r(2A + i)nr(M, + i)(-i)^ 

i^i 

r(A + |)r(|M|+Ar + 2A + l) 



{B.20) 



X < 



X 



E 



2l/2(iV-|fc|) 



h 

1 "'p 



T{\M\ + l{N-\k\ + \ + ^) 



N = \k\ (mod 2) 

^) ( -M, + ki,...,-Mp + kp f l-Uj-n 
i -|M|-i(7V-|^|), + l 1 2 



p 

X E n i^J ■ n^., )S (n^. \ U^=i^fc.) ' ■ 

n Afc. = 0, 

p p 

Here we have used the notation (B.15), (143) and \M\ — ^ Mj, \k\ = ^ kj. 

3=1 j=i 

Prom formula (B.20) we obtain an expression for the multipoint correlation 
functions of the classical Heisenberg 0(m)-model on a complete graph (and 
hence on an arbitrary graph). The answer is a relation of generalized hy- 
pergeometric series which is too cumbersome and requires further study and 
simplification. 

In conclusion we note that identity (B.18) gives an alternative means of 
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computing the statistical sums {ip \ iIj)q and (ip \ ip)rcg (see (133) and (135)) of 
the classical Heisenberg 0(m)-model on open and closed chains respectively. 
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